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Abstract 

Recently, C. Imbert & R. Monneau study the homogenization of coercive Hamilton- 
Jacobi Equations with a w/e-dependence : this unusual dependence leads to a non-standard 
cell problem and, in order to solve it, they introduce new ideas to obtain the estimates on the 
oscillations of the solutions. In this article, we use their ideas to provide new homogenization 
results for "standard" Hamilton-Jacobi Equations (i.e. without a u/e-dependence) but in 
the case of non-coercive Hamiltonians. As a by-product, we obtain a simpler and more 
natural proof of the results of C. Imbert & R. Monneau, but under slightly more restrictive 
assumptions on the Hamiltonians. 



Introduction 

In this article, we are interested in homogenization problems for first-order Hamilton- 
Jacobi Equations. The originality of this work is to provide results in the case of non- 
coercive Hamiltonians and applications to non-standard problems. Before describing more 
specifically our contributions, we want to point out that most of the new ideas used in 
this paper are borrowed from Imbert & Monneau [20] who study the homogenization of 
(coercive) Hamilton-Jacobi Equations with a M e /e-dependence, namely 

u\ + Hie^x, e~ x u E , Du E ) = in R n x (0, +oo), (1) 

where if is a continuous and coercive Hamiltonian. In fact, the starting point of the 
present work and one of its main motivations was to solve such problems, maybe under 
more restrictive assumptions than in [20] but (i) with simpler and more natural methods 
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and (ii) with a clearer explanation of the involved phenomenas. We will come back later 
on (pQ) and, in particular, we will explain the connections between the standard problems 
we are first dealing with and these non-standard problems. 

In order to emphasize the main new ideas in our approach, we are not going to try 
to consider the most general framework but we restrict ourselves to a model case which 
carries the main difficulties. More precisely, we are going to study the limit as e — > of 
the solution U £ (x,y,t) of the Hamilton- Jacobi Equation 

U £ + F(e~ 1 x, e~ l y) e~H, D X U £ , D y U £ ) =0 in R rt+1 x (0, +oo), (2) 

where x G R n , y G R, t G (0, +oo), F(x,y,t,p x ,p y ) is a continuous function on R 2n+3 
which is Z n -periodic in x and 1-periodic in y and t. More precise conditions on F will be 
given later on but we point out that F is assumed to be coercive with respect to p x but 
not with respect to p y and this is the key new point. In order to give a flavour of the type 
of Hamiltonians we are able to consider, a typical example is 

F(x,y,t,p x ,p y ) = a(x,t)\p x f + b(x,y,t)\p y \ - f(x,t) , 

where a, b, f are Lipschitz continuous function with the right periodicity in x, y, t with 
a(x, t) > r] > in R n x [0, +oo) and j3 > 1. In particular, no sign condition is imposed on 
b, it may change sign and even the term b(x,y,t)\p y \ can be replaced as well by a linear 
term b(x, y, t)p y . 

In order to understand the contribution of this work, we recall that the theory of 
homogenization for first-order Hamilton- Jacobi Equation started with the famous unpub- 
lished work of Lions, Papanicolaou & Varadhan [27] who completely solve the problem 
in the case of time- independent, periodic and coercive Hamiltonians i.e. for the above 
equation when 

F(x,y,p x ,p y ) — > +oo when \p x \ + \p y \ — > +oo, uniformly w.r.t x and y. 

Then, to the best of our knowledge, such an assumption was used (one way or the 
other) in all the works concerning the homogenization of Hamilton- Jacobi Equations : 
for more general periodic situations (cf. Ishii|23j). for problems set in bounded domains 
(cf. Alvarez [T], Horie & Ishii [22] ) , for equation with different structure (cf. Alvarez & 
IshiijS]), for deterministic control problems in L°° (cf. Alvarez & Barron [5]), for almost 
periodic Hamiltonians (cf. Ishii[24j) and for Hamiltonians with stochastic dependence (cf. 
Souganidis[28]). 

The reason why this coercivity assumption plays a central role, is that it allows to 
solve the so-called "cell problem" which provides the homogenized equation, i.e. the 
homogenized Hamiltonian. When F is independent of time, this problem, which turns 
out to be an ergodic problem, consists in solving the pde 

F(x,y,D x v+p x ,DyV+ Py ) = X in R" +1 , 

and in showing that, for any (p x ,p y ), there exists a unique constant A = F(p x ,p y ) such 
that this equation has a bounded (periodic) solution. In general, solving this cell problem 
is the main difficulty and, in most cases, the fact that F is, indeed, the homogenized 
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Hamiltonian, follows rather easily. It is worth pointing out, anyway, that, in order to 
have a comparison result for the limiting equation, one has to prove that F satisfies 
suitable properties and this may be a difficulty in some cases. We refer the reader to 
Section H] for a discussion in this direction. 

In this framework, the role of the coercivity is rather clear since it provides an easy 
bound on Dv = (D x v, D y v) once one knows that F(x, y, D x v + p x , D y v + p y ) is bounded 
and this is a key argument to solve the cell problem. Up to now (again, to the best of our 
knowledge), no work (except perhaps partially |28j) succeeded to find an argument to solve 
this cell problem by pde methods in a general way which avoids the coercivity assumption 
on F (and actually one of the main question is to find the right set of assumptions to do 
it). Furthermore, when F depends on t, one has to find a space-time periodic solution of 

Vt + F(x, yj t,D x v + p Xj D y v + Py) = F(p x ,p y ) inl n+1 xl, (3) 

and the assumptions to do it seem even stronger (cf. Souganidis and the author [T7]). 

However it is worth pointing out that the pure pde approach we describe above for 
both homogenization and ergodic problems can be replaced, in some cases, by an approach 
using either optimal control (if F is convex) or differential games methods : in this 
framework, results do exist for non-coercive Hamilton- Jacobi Equations, the coercivity 
assumption being replaced by either controllability or non-resonance assumptions. We 
refer the reader to Alvarez and Bardi [21 El H] , Arisawa [71 [5] , Arisawa and Lions [H] and 
Artstein and Gaitsgory [10J, Bardi [Uj for results in this direction. 

In order to solve problems like ([3|) in a general way and once for all {p x ,p y ), we consider 
the ergodic problem : find a constant A such that the equation 

w t + G(x, y, t, D x w, D y w) = A inl n+1 xR, (4) 

has a bounded (periodic) solution. In fact, with our approach, we obtain a weaker but 
sufficient result which can take two possible forms : we can prove that there exists a 
unique constant A such that we have either approximate continuous (but not necessarily 
Lipschitz continuous) periodic solutions or exact, possibly discontinuous, periodic sub and 
supersolutions. This result is the first main contribution of the paper and this is where we 
use in a key way the ideas of Imbert & Monneau [20J ; an unusual feature of the proof is 
the estimate of the oscillation max r +i X |(«)) — min R n+i xR (w) which replaces the classical 
gradient estimate and relies on two very original arguments. 

Once (HI) is solved, the result for the homogenization problem follows by using the 
usual arguments : indeed, it is enough to (essentially) apply the result to for any 
(Px,Py), then to obtain suitable properties for F and finally to use the "perturbed test 
function's method" of Evans [HI [19], even if, to prove the convergence, we have to intro- 
duce an additional argument to take care of the rather weak properties we impose on F 
in the variables x, t and p x . We refer the reader to the papers of Alvarez & Bardi [21 [3] for 
a general approach of singular perturbation problems for fully nonlinear, possibly degen- 
erate, elliptic or parabolic equations; this approach applies to homogenization problems, 
including in the case of first-order Hamilton- Jacobi Equations, and, in particular, the au- 
thors clarify the connections between ergodic and homogenization problems in a general 
setting. 
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Now we come back to (pQ) which seems to have nothing to do with the above (rather) 
classical problem. Roughly speaking, the connection consists in considering as the 
equation of the motion of the graph y = u e (x, t) and to introduce the associated level set 
equation which is an equation in R n+1 . Since, intuitively, the level set equation should 
play a role only on the moving hypersurface, this leads to replace e u e (x, t) by s~ x y 
and to transform (CQ) in a standard (but a priori non-coercive) homogenization problem. 
The rigourous justification of this formal (but convincing) argument is easy by following 
the article of Giga & Sato [21] which is exactly describing the way to do it for first- 
order equations (see also the paper of Biton, Ley & the author[T3j for motion by Mean 
Curvature). Compared to the article of Imbert & Monneau [20J, this approach leads to 
more restrictive assumptions on H but to a far simpler proof and to more natural and 
understandable arguments. Typically the results of [20] apply for all H of the form 

H{x,y,t,p) = a{x,t){\p\ 2 + l)^ 2 + g{x,y,t) , 

where a, g are continuous periodic functions with a > in R n x R and < (3 < 1, while 
here we have to consider the case (3 = 1. Finally, it is worth pointing out that Imbert & 
Monneau [20] use also an extension to R ra+1 by adding a extra variable but with no clear 
interpretation of this new variable which is just a trick in the proof. 

The paper is organized as follows : in the first section, we study the ergodic problem 
01]); this is the occasion to present the key assumptions on G (and therefore on F) which 
describe the structure of the problem allowing the lack of coercivity. In Section [2] we 
apply this result to the cell problem and we solve the homogenization problem; to do 
so, we show that F(p x ,p y ) is a continuous function, which depends also continuously on 
F. Then, in Section [3] we consider the problem (TTJ: we describe the Giga & Sato [2"Tj 
approach and the assumptions on H to apply it; the conclusion then follows from the 
results of Section [2] Finally, we provide remarks on the possible variants and extensions 
of the above results in Section HI 

Acknowledgments : The author wishes to thank Martino Bardi for interesting com- 
ments on the first version of this paper and, in particular, for pointing out some missing 
references. He also thanks the anonymous referees for their careful reading of the paper 
and their very constructive remarks. 

1 The non-coercive ergodic problem 

In order to state the result, we have to impose two kinds of assumptions: first some 
basic assumptions which ensure existence, uniqueness and comparison properties for the 
problems we are going to introduce and then structure conditions which allow to obtain 
the estimates we need to solve the ergodic problem. Again, to focus on the main ideas 
of the proof, we are not going to provide the most sophisticated assumptions but the 
simplest relevant ones. 

(HI) G(x,y,t,p x ,p y ) is a continuous function in R 2n+3 , Z n -periodic in x, 1-periodic in y 
and t. 
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(H2) Either G is independent of t and G(x,y,p x ,p y ) — > +00 when \p x \ — > +00, uniformly 
w.r.t. (x,y) G R n x R and jpj < i? for any R > 0, or G depends on i and there exists 
constants C\, C2 > such that 

G(x, y, t, p X} 0) > d \p x \ - C 2 in f x R x 1 x R rt . 

(H3) The function G is locally Lipschitz continuous w.r.t. y, t, p y and there exists / G R 
and C3, 6*4,(75 > such that, for almost every £ = (x,y,t,p x ,p y ) G R 2n+3 

IAM0I < Cab, + J| , |AG(0I < C 4 (l + b,| + |G(0I) , < Cs- 

Our result is the following 
Theorem 1.1 Under assumptions (H1)-(H3), we have 

(i) For any a > 0, there exists a unique continuous, space-time periodic solution w a of 

< + G(x, y, t, D x w a , D y w a ) + aw a = in R" +1 x R , 
which is independent oft if G is independent oft. 

(ii) There exists a constant K = K{G) depending only on G through C := | \G(x, y, t, 0, 0)| | 
I and the coercivity assumption (H2), such that 

max w a — min w a < K . 

In particular, if G depends on t, then K = K(Cq, C\, C2, 1). 
(Hi) There exists a unique constant A such that the equation 

w t + G(x,y,t,D x w,D y w) = \ in R n+1 x R , (5) 

admits bounded, possibly discontinuous, sub and super solutions . 

(iv) If Vq is bounded uniformly continuous in R n+1 ; there exists a unique solution v of the 
initial value problem 

v t + G(x, y, t, D x v, D y v) = in R n+1 x (0, +00) , (6) 

v{x,y,0) = v {x,y) ml n+1 , (7) 

which is bounded and uniformly continuous in R n+1 x [0, T] for all T > 0; moreover we 
have 

lim t _1 t>(x, y, t) = — A , 

the limit being uniform in R n+1 . 

We have chosen to state Theorem 11.11 in that way for several reasons : first, even 
if the result (i) seems to be classical, it is not the case because the assumptions (Hl)- 
(H3) are not completely standard and it may not be so clear that the iy a -equation has a 
unique solution. Indeed, we only assume G to be continuous in x, which is a very weak 
assumption, but this is compensated by the coercivity of G in p x . To prove (i) which 
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is the first (slight) difficulty of the proof, we borrow arguments from P.L. Lions and the 
author [J5] and the author [13J. Then (ii) is the main step but (curiously) not the most 
technical one. (hi) is one of the possible conclusion but if one insists to use continuous 
"correctors", then (ii) implies that, for any 5 > 0, if a is small enough, we have 

A - 5 < < + G(x, y, t, D x w a , D y w a ) <X + 5 in R n+1 x R , 

since, as usual and by using (ii), A is the uniform limit of — aw a . Finally (iv) is an other 
classical characterization of the ergodic constant A. 

Proof : We provide the proof in the case when G actually depends on t, the other case 
being simpler (see the remark at the end of the section). 

1. We start by (i). We just sketch the proof and refer to [15], [13] and [12] for the (easy) 
details. As it is classical in viscosity solutions' theory, the existence of w a relies on having 
a Strong Comparison Result for the equation, the existence being obtained by the Perron's 
method extended to the framework of viscosity solutions by Ishii [22]. Here, because of 
the periodicity of G, Perron's method can be applied in the set of periodic subsolutions 
(by remarking that if w is as a subsolution, w(x,y,t) := sup{w(x + l,y + k,t + m, I G 
Z n , k, m G Z} is still a subsolution), we can prove the comparison result only for periodic 
sub and supersolution which avoids problems with the unboundedness of the domain. 

To prove this comparison result, we argue in the following way: if u is an use pe- 
riodic subsolution and v a lsc periodic supersolution of the equation, we approximate 
maxp+2 (u — v) by max R 2n+4 (\1/) where \I/ = ^ v ,/3 is given by 

W(xi, x 2 , yi, 2/2, t, s) := u(xi, y h t) - v{x 2 , y 2 , s) — 2 — 2 — 2 — , 

where rj, [3 > are small constants devoted to tend to 0. 

Because of the periodicity of u and v, the use function \1/ achieves its maximum at 
some point, which is still denoted (xi, x 2 , yi,V2, t, s) for the sake of simplicity of notations 
(but which depends on rj and (3). We set 

2{x 1 -x 2 ) , 2( yi -y 2 ) ^ 2(t-s) 



Px ■= 2 , Py ■= ~T 2 , Pt 



r] 2 ' " y (3 2 ^ (3 2 

The viscosity solutions inequalities for u and v read 

p t + G(x 1 ,y 1 ,t,p x ,p y ) +au(x 1 ,y 1 ,t) < , (8) 

Pt + G(x 2 , 2/2, s, p x , p y ) + av(x 2 , 2/ 2 , s) > . 
At this point we recall the well-known fact that 

(We point out here that this is true since we deal with time-space periodic sub and 
supersolutions). 
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This implies, in particular, that p t = o{(3~ 1 ) and since u and v are bounded, we may 
as well assume that, for f3 small enough, the above viscosity inequalities holds with G 13 = 
min(/3 _1 , max(G, — which has the same properties as G except that \D t G^(^)\ < 

c^i + bj + r 1 ). 

Subtracting the above viscosity sub and supersolutions inequalities and using (H3) 
then yield 

a(u(xi,yx,t) - v(x 2 ,y2,s)) < G (x 2 ,yi,t,p x ,p y ) - G p {x 1 ,y 1 ,t,p x ,p y ) 

+C 3 \p y + l\\yi - y 2 \ + C 4 \t - s|(l + \p y \ + (3- 1 ) . 

The argument is then the following : we first fix (3 and let r\ tend to 0. Because of (H2) 
and inequality (jEJ), the ^-gradient termpa, remains bounded since pt,p y remains bounded ; 
furthermore, since u and v are periodic, we may also assume that (x\,X2,y\,y2,t,s) re- 
mains in a compact subset of IR 2n+4 ; therefore the passage to the limit 77 — > can be done 
by using only the continuity of G@ (or, equivalently, of G), all terms being convergent up 
to a subsequence. 

Therefore, we are left with 

a(u(x 1 ,y 1 ,t) -v(x 2 ,y 2 ,s)) < C 3 \p y + l\\y x - y 2 \ + C 4 \t - s|(l + \p y \ , 

and the conclusion follows by letting (3 tend to : the terms of the right-hand side converge 
to because of (jUj), while we know that u(xi, y±, t) — v(x 2 , y 2 , s) — > max R n+2 (u — v); we 
reach the conclusion that max K n+2 (u — v) < and the proof of (i) is complete. 

2. We further remark that, either by construction or by comparison with constants, we 
have 

min [-G(x, y, t, 0, 0)] < aw a < max [-G(x, y, t, 0, 0)] . (10) 

x,y,t x,y,t 

In particular, ||aw Q ||oo < Co- 

3. Now we turn to (ii). The idea of Imbert & Monneau consists in estimating separately 
the behavior of w a in (x, t) and in y. For the estimate in (x, t), we introduce the function 

w(x, t) := max w a (x, y, t) , 



which, using (H2) and the estimate on aw a given in Step 2, is a space-time periodic 
subsolution of 

w t + C 1 \D x w\-{C + C2) <0 inl n xl. 
But the Oleinik-Lax Formula shows that, for any x, t and for any s < t 

w(x, t) < min [w(z, s)] + (C + C 2 )(t - s) , 

\z—x\<Ci(t—s) 

and therefore, for any S > 

w(x, t) < min <^ min [w(z, s)} + (C + C 2 ){t - s) 

0<t—s<S \\z—x\<Ci(t—s) 
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Choosing (x, t) such that w(x, t) = max PxK w and choosing S large enough in order that 
the set {(z, s); < t — s < S , \z — x\ < Ci(t — s)} contains a whole period, we deduce 
that 

max w — min w < K\ , (11) 
where Ki depends only on C , C\ and C 2 . 

4. Next we consider the behavior of w a in y. We are going to show that D y (w a (x, y, t)+ly) 
has the same sign as I. We only do it in the case / < 0, the case / > being treated in 
an analogous way. If I = then w a is independent of y : indeed, since G(x,y,t,p x ,0) is 
independent of y by (H3), by using the same arguments as in the proof of (i), one can 
build a solution for this new Hamiltonian which depends only on x and t and which is 
also a solution of the G equation. Therefore, by uniqueness, w a = w a (x,t). 
To prove that y h- > w a (x, y, t) + ly is nonincreasing for / < 0, we consider 

M := max{w a (x, yi ,t) - w a (x,y 2 ,t) + l( Vl - y 2 ); x E R n , t E R, yi > y 2 } ■ 

We have to show that this maximum is nonpositive. To do so, we argue by contradiction 

assuming that M > and we introduce the function 

/ . n a , a , \ (^i-x 2 ) 2 [{yi-y 2 y} 2 
X{x 1 ,x 2 ,y 1 ,y 2 ,t,s) := w (x 1 ,y 1: t)-w {x 2 ,y 2 ,s) — — 

The function x nas indeed a maximum since w a is continuous and periodic (hence 
bounded) and since / is strictly negative, the term l{yi — y 2 ) controls (yi — y 2 ) + which 
therefore remains bounded. 

Denoting by (x±, x 2 ,y\, y 2 , t, s) a maximum point of x an d 

„_2(t-s) _ 2{ Xl -x 2 ) _ 2{ yi -y 2 )- 
i Vx o ' Py 



p 2 > r* ^ ry - p 2 

we have the viscosity inequalities 

a + G(x 1 ,y 1 ,t,p x ,p y ) + aw a (x 1 ,y 1 ,t) < , 
a + G(x 2 ,y 2 ,s,p x ,p y ) + aw a (x 2 ,y 2 , s) > . 

Using (H3), we have 

\{yi-y2)-? 



\G(x 1 ,y 1 ,t,p x ,p y ) -G(x 1 ,y 2 ,t,p x ,p y )\ < 2C 3 - 



and subtracting the above inequalities, we are left with 

a(w a (x 1 , y h t) - w a (x 2 , y 2 , s) + Z(?/i - y 2 )) < 

G(x 2 ,y 2 ,s,p x ,p y ) - G(x 1 ,y 2 ,t,p x ,p y ) + 2C 3 — ^ - +al(y 1 -y 2 ) . 

To conclude, we argue as in the proof of the comparison result : first, by periodicity we 
may assume that (xi, y±, t) and (x 2 , y 2 , s) remain in a compact subset of M n xRxR; next, 
by classical arguments, all the penalization terms tend to 0. In particular, (yi—y 2 )~ = o{(5) 



8 



and this yields a.l{y\ — y<i) < o(/3) since I < 0. Therefore we first fix j3 and let r\ tend to 
and then we let (5 tend to 0, using the same arguments as above. This gives the inequality 
aM < and the conclusion. 

5. To prove (ii), it is enough to put together the two above informations and again we 
do it only for / < 0: we consider two points (x,y,t) and (x',y',t') such that w a (x,y,t) = 
maxRn+2 w a and w a (x', y', t') = minjRn+2 w a , and we denote by y the smallest point z > y' 
such that w(x',t') = w a (x' , z,t'); we recall that w a is periodic and therefore y — y' < 1. 
We have 

w a (x,y,t) 





w(x, t) 




< 


w(x , t') 




< 


w a (x',y, 


t') + Ki 


< 


w a (x',y> 


,t')+l(y'-y)+K 


< 


w a (x',y' 


,t')-l + K 1 . 



which concludes the proof of (ii). 

6. For the proof of (iii), we notice that w a — min w a is bounded and so is a min w a since 
||a!iu a ||oo < Co- Without loss of generality, we may assume that —a min w a converges to 
a constant A and, along the same subsequence, the half-relaxed limits of w a — min w a , 
denoted by w and w, provide respectively a bounded subsolution and a bounded super- 
solution of ([5]). 

Now we show that there exists a unique constant A such that ([5]) has a bounded 
subsolution and a bounded supersolution. First, we assume that there exists A such that 
§5§ has a bounded subsolution u. By comparison result for the initial value problem (a 
result proved in a similar way as in the stationary case), using that u — Xt and w — Xt are 
respectively sub and supersolution of equation ([6]), we deduce 

max (u — w) (x, y, t) < max (u — w) (x, w, 0) + t(X — X) . 

Since u and w are bounded, by dividing by t > and letting t tend to +oo, we obtain 
A — A > 0. The reverse inequality is obtained by the same argument assuming that there 
exists a bounded supersolution for the same constant. 

7. Finally we prove (iv). The proof of the existence and uniqueness of the solution 
v follows along the line of the proof of (i) with few additional classical arguments and 
therefore we skip it. In order to prove the convergence property, we use w and w defined 
in step 6 above. Since vq, w and w are bounded, there exist constants c, c such that 



c + w(x, y, 0) < vo(x, y) < c + w(x, y, 0) in 



D n+1 



But c + w(x, y, t) — Xt and c + w(x, y, t) — At are respectively sub and supersolution of 
and, by comparison, we deduce 



c + w(x, y, t) - Xt < v(x, y,t) < c + w(x, y, t) - Xt in R" + x (0, +oo) . 

Dividing by t and using the fact that W and w are bounded provide the property we 
wanted to prove. And the proof of Theorem 11.11 is complete. 
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Remark 1.1 In the case where G is independent oft, w a is also independent oft since 
w a (x, y, t) andw a (x, y, t+h) are both solutions of the equation for any h G M and therefore, 
by uniqueness, w a (x,y,t) = w a (x,y,t + h) for any h G R. But maybe a more simple 
approach in this case is to solve the stationary equation directly and not to consider the 
evolution equation. In the above proof, the only change concerns the estimate obtained 
through W = w(x) : indeed W is a subsolution of 

min G(x, y, Dw, 0) < C in R n , 

and this inequality together with the coercivity asssumption on G, shows that \ \Dw\ |oo < K, 
where K depends only on G through Cq and the coercivity. The function w being Z n - 
periodic in x, this gives the bound on the oscillation and even a more precise information. 

2 Homogenization of the non-coercive equation 

In order to state and prove the homogenization result, we first need to identify the ho- 
mogenized Hamiltonian and its properties. 

Lemma 2.1 Assume that F satisfies the assumptions (H1)-(H3) with I = in (H3). 
For any (p x ,Py) G M n x M, there exists a unique constant F(p x ,p y ) such that ([5|) has a 
bounded space-time periodic solution. The function F(p x ,p y ) is a continuous function of 
(Px,P y )- Moreover, if (F k ) k is a sequence of functions satisfying the same assumptions as 
F and uniformly in k, and which converges locally uniformly to F on IR 2n+3 ; then F k — > F 
locally uniformly in M n x M. 

Proof : The first part of the result is an immediate consequence of Theorem 11.11 applied 

to 

G(x, y, t, q x , q y ) = F(x, y, t, q x + p x , q y + p y ) . 

The continuity of F comes from all the estimates we have : first we have a bound on 
F(p x ,p y ) which is coming from (fTUj) (recall that —aw a converges uniformly to A), namely 

min[F(x,y,t,p x ,p y )\ < F(p x ,p y ) < max[F(x,y,t } p x ,p y )} . 

x,y,t x,y,t 

Then we have also the bounds on the oscillations of the associated sub and supersolution 
v, v, typically through K\ given by f TTTT) in the time dependent case, which depend only 
on {p x ,p y ), F and I = p y . 

Indeed, if the sequence ((p x ,P y ))k is converging to (p x ,P y ), then, up to a subsequence, 
F(p x ,p y ) converges to /x G R. But if v k and v k are bounded sub and supersolution associ- 
ated respectively to (p x ,p y ), then they are uniformly bounded and the half-relaxed limit 
method provides v and v which are the sub and supersolution associated to (p x ,p y ) and 
the constant /i. By the uniqueness result of Theorem 11.11 (iii), this implies \i = F(p x ,p y ) 
and since this is true for any converging subsequence of (F(p x ,p y ))k, the continuity is 
proved. 

Finally the proof for the F relies on the same type of arguments and therefore we 
skip it. 

The homogenization result is then the 
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Theorem 2.1 Assume that F satisfies the assumptions (H1)-(H3) with I = in (H3). 
If the U £ are continuous viscosity solutions of ^) associated with the initial data U £ (x, y, 0) = 
Uo(x,y) in M. n x R where Uq is bounded uniformly continuous in M. n+1 , then, as e — ► 0, 
U £ — > U locally uniformly in M n+1 x (0, +oo), where U is the unique solution of 

U t + F(D X U, D y U) =0 in « n+1 x (0, +oo), (12) 

with the same initial datum U 



Proof : We just sketch it since it follows essentially by classical arguments. We first 
remark that, if A : = max(||£/o||ooj | \F(x, y, t, 0, 0)| |oo), then —A(l + t) and A(l + t) are 
respectively sub and supersolution of the problem and, by comparison, we have 

-A(l +t)< U £ (x, y, t) < A(l + t) in R n+1 x (0, +oo) ; 

therefore, the U £, s are uniformly bounded. 

Then we use the half-relaxed limit method and set u := limsup* U £ , u = liminf* U £ . 
We show how to prove that u is a subsolution of (112D . the proof that u is a supersolution 
being analogous. 

If cf) is a smooth function and (x, y, t) is a strict maximum point of u — <p, we use the 
perturbed test-function's method in the following way: we set p x = D x cj)(x,y,t), p y = 
D y <f)(x,y,t) and let V be a bounded space-time periodic supersolution of ([3]) associated 
to (p x ,Py) and to the function 

F k (x } y } t } q X} q y ) = min {F(x', y, t r , q' x , q y ) ; \x' - x\, \t' - t\, \q' x - q x \ < &T 1 } , 

for k > 1. 

We introduce the function 

-i / -i / — \ x ~ x I 1 2/ 2/ I 1^ — t I 



U £ (x, y, t) - cj)(x, y, t) - eV(e~ x , e" V, e"Y) - 



S 2 S 2 S 2 



It is easy to see (and standard) that, as e — > and 5^0 with 5 e, the maximum points 
of this function converge to (x,y,t,x,y,t). Performing standard computations (which we 
recall in the Appendix for the reader's convenience) with fixed k leads to the inequality 

<f>t{x,y,t) +F\D x <t>(x,y,T),D y( t>(x,y,t)) < , (13) 

and we use Lemma [27X1 to conclude by letting k — > oo. 

We just point out that the F^-trick allows to take into account the rather weak condi- 
tions we impose on F(x, y, t, q x , q y ) in x, t and q x : indeed, if (x, y, t, x', y', t') is a maximum 
point of the above function, we have, on one hand, \D x <f)(x, y, t) — D x (fi(x, y,t)\ < k~ l for 
e and 5 small enough because (x, y, t) — > (x, y, t) when e, 5 — >• 0; on the other hand, since 
the terms ^| ■ , ^~ s l ■ are bounded, then \e~ l x' — e~ l x\, \e~ l t' — e~ l t\ < k^ 1 if 5 <^ e is 
small enough. Therefore 

F k (e' 1 x', e~ x y\ e~H', D x <p(x, y, t) + q x , q y ) < F{e' 1 x, e^y', e~H, D x <p(x, y, t) + q x , q y ) , 
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which is the inequality needed in the proof because the dependence in y, y' and q y is taken 
into account by (H3). Without this argument, we would face a difference of F-terms 
which is not a priori small. A more detailed argument is given in the appendix. 

We end the proof by remarking that we have a Strong Comparison Result for the F 
equation since F is continuous and just depends on (p x ,Py) and this allows to compare 
the subsolution u and the supersolution u to get the complete answer (see, for example, 
Ley[2S])- This point may be a difficulty in more general cases and we refer the reader to 
Section H] for remarks in this direction. 

3 Homogenization of the u e je - equation 

In this section, we consider (pQ) and even a more general time dependent one, namely 

u E t +H{e~ 1 x,e~ 1 u £ 1 e~H,Du e )=Q in R™ x (0, +00), (14) 
together with the initial data 

u £ (x,0) = u (x) inR n . (15) 

where uq is bounded, uniformly continuous in R n and the Hamiltonians H e C(R n xlx 
R n ) satisfies the assumptions 

(H4) H(x,u,t,p) is a locally Lipschitz function in R n x R x R x R 71 , Z n -periodic in x, 
1 periodic in u and t. Moreover there exists a constant C such that, for almost every 
(x,u,t,p) in R" x R x R x R" 

\D t H(x,u,t,p)\ < C(l + \H(x,u,t,p)\) , \D u H(x,u,t,p)\ , \D p H(x,u,t,p)\ < C. 

(H5) H is coercive in p, i.e. 

H(x,u,t,p) — > +00 as \p\ — > +00 uniformly w.r.t. x G R n , u, t £ R. 

In addition to these assumptions, we have to add a "geometrical" assumption which 
ensures that (fT4l) is the graph-equation associated to a level-set equation in R n+1 . A priori 
this requires only 

(H6w) There exists an Hamiltonian Hoc satisfying (H4)-(H5) such that 

iJoo(x, m, t, p) = lim sH(x, it, t, s^ x p) . 

but we use below a slightly more restrictive one 

(H6s) There exists a constant C such that \D p H(x,u,t,p) ■ p — H(x,u,t,p)\ < C a.e. in 
R x R n x R™. 

"w" is for weak and "s" is for strong. Below we just use the notation (H6) for (H6s) 
but keeping in mind the notation introduced in (H6w). 

Before providing the result, we point out that both (H6w) and (H6s) hold for 
H(x, t, u,p) = c(x, t)\p\ +g(u, t) if c, g are Lipschitz, periodic functions with c(x, t) > r\ > 
in R n . Notice that g may change sign or be negative. 

The result is the following 
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Theorem 3.1 Assume that H satisfies (H4)-(H6). Then the sequence {u £ ) £ converges 
locally uniformly to a function u which is the unique solution of an equation of the type 

u t + H(Du) = in W 1 x (0,+oo), (16) 

with 

u(x,0)=u (x) mW. (17) 

The Hamiltonian H can be interpreted as an homogenized Hamiltonian for a standard 
homogenization problem in dimension n + 1 and is characterized by the following property: 
A = —H{p) is the unique constant such that, if w is the solution of 

w t + H(x,w,t,Dw) = mR n x(0,+oo), (18) 
w(x, 0) = p ■ x + w (x) inR n , (19) 

where Wq is a bounded uniformly continuous function, then w(x,t) —p- x — Xt is uniformly 
bounded in M. n x (0, +oo). In particular 

—H(p)= lim t~ 1 w{x,t) , 

the limit being uniform in M. n . 



Proof : 

1. From M. n to M ra+1 : we interpret the equation for u e as the equation of the motion 
of a graph and we introduce the level-sets equation in M n+1 . To do so, we use (a slight 
modification of) the approach used by Giga & Sato [21] for Hamilton- Jacobi Equations 
(see also Biton, Ley and the author [H] for the Mean Curvature Equation). The function 

U e (x, y, t) := u £ (x, t) — y solves 

U £ t +\D y U £ \H(e- l x,e- l U £ + e- l y,e~H,\D y U £ \- l D x U £ ) = in R n+1 x (0,+oo), 

U £ (x,y,0) = u Q (x) -y in R n+1 . 

In the sequel, it is convenient to use the notation X — (x, y), P — (p x ,Py) an d to 
introduce the Hamiltonian F defined by 

F(X,t,P) :-- 



\Py\H(x, y, t, \p y \ l p x ) if p y ^ 0, 
(x,y,t,p x ) otherwise. 



By the assumptions (H4)-(H6) on H, tedious but straightforward computations show 
that F satisfies (H1)-(H3) ((H6s) provides, in particular, the Lipschitz continuity inp^). 
It is worth pointing out that, in (H2), since F is homogeneous of degree 1, the coercivity 
of F in p x or the stronger property required in the t-dependent case are, in fact, equivalent. 

We finally consider x e (X, t) := l^jje^.^> y where, here and below, 1^ denotes the 
indicator function of the set A. Following [16j, it is easy to see that it solves 

X £ t + F{e- l X, e~H, D X £ ) = in R n+1 x (0, +oo), 
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X e (x,0) = l {uo(:E) _,>o } inR n+1 . 

Indeed the result is obvious if either U £ (x,t) > or U £ (x,t) < because x e is constant 
in a neighborhood of (x,t) and the equation is a trivial consequence of the fact that 
F(X,t,0) = for any X and t. Hence we are left with the case when U e (x,t) = 0, 
where the arguments of [16] provide the answer by using the continuity of U e . In fact, in 
the equation, the dependence in U e is treated as a dependence in X and t ; the change 
of variables of [16] is done only at the level of the derivatives. This explains why the 
dependence in U e disappears from this equation. 

2. Homogenization in R n+1 . The x e -equation can be homogenized by using almost readily 
Section [2j the only slight difficulty being the discontinuous initial data and solutions. In 
fact, if x '■= limsup*x e , X '■= liminf*x e , u := limsup*u £ , u := liminf*^, we have the 
following: x and x are respectively viscosity sub and super solutions of 

Xt + F(D X ) = inR n+1 x (0,+oo), (20) 

with 

{u (x)- y >o} m K , 

and 

X(X,t) = l {lI ( x ,t)-y>0} , X(X,t) = l { u(x,t)-y>0} IB X (0, +Oo). 

We also point out that F is homogeneous of degree 1. 

3. Reconstruction of the result in R n . Let u : M n x (0, +oo) —>■ R be the solution of 

Ut + F(D x u, -1) = in R n x (0,+cx)), 

u(x, 0) = u (x) inR". 

By the arguments of Soner, Souganidis and the author [16] . since u(x, t) — y is the solution 
of the level set equation ( l20i) associated to the initial data uo(x) — y, we have 

X_(X,t) > l{ u ( x ,t)-y>o} and x(X,t) < l{«(*,t)- w >o} in x (0, +oo). 

Indeed, it is proved in [16] that l{ u (a;,t)-y>o} an d ^{u(x,t)-y>o} are respectively the minimal 
supersolution and the maximal subsolution associated to the initial data H-{ Uo (x)-y>o}- By 
the form of x and x given in Step 2, this yields 

u(x, t) < u(x, t) < u(x, t) < u(x, t) in R n x (0, +oo). 

Therefore u = u = u in R n x (0, +oo) and this proves the first part of the result with 
//(/,): r(/,-l). 

4. To prove the last part of the result, we apply the same approach : if U : R n+1 x 
(0, +oo) — » R is the solution of the level set equation 

U t + F(X,t,DU) = in R n+1 x (0,+oo), (21) 

with initial datap-x+w (x)— we first have {(x, y, t); w(x,t) = y} C {(x,y,t); U(x,y,t) = 
0}. 
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On the other hand, if V\ and v 2 are bounded space-time periodic functions which are 
respectively an use subsolution and a lsc supersolution of ([3]) associated to P = (p x , — 1), 
then the functions v 1 (X,t) + P ■ X - F(P)t, v 2 (X,t) + P ■ X - F(P)t are respectively 
sub and supersolution of fl2TT) ; since, we may assume in addition that max R n+2 v\(X,t) < 
wo < min R n+2 v 2 (X, t) in M n+1 , we have by using a comparison result 

Vi(X,t)+P-X- F(P)t < U(X, t) < v 2 (X, t) + P ■ X - F(P)t in R" +1 x (0, +oo) . 

Notice that, since H satisfies (H4) and (H6), F is Lipschitz continuous in P, uniformly 
in X and we can use a property of "finite speed of propagation" for the level set equation, 
which simplifies the comparison arguments and avoids problem with the unboundedness 
of the domain. 

If w(x, t) = y, then U(x,t) =0, P ■ X = p ■ x — w(x, t) and this inequality gives 
Vi(X,t)+p-x-w(x,t)-F(P)t < < v 2 (X,t)+p-x-w(x,t)-F(P)t in R n x (0, +oo) , 
or equivalently 

Vi(X,t) < w(x,t) -p-x + F{P)t < v 2 (X,t) inW 1 x (0, +oo) . 

Since H(p) = F(p,—1) and v±, v 2 are bounded, this gives the first part of the result. 
Then, dividing by t and letting t — > oo, we obtain the second one. And the proof is 
complete. 

4 Extensions and Remarks 

The results of the preceding sections can be extended in several directions. 

First, in Section [2], we can consider stationary homogenization problems. In fact, as 
we suggest it in the introduction, these problems are simpler provided that we have a 
suitable dependence in u e . 

Next, the one-dimensional variable y can be replaced by a multi- dimensional variable 
(yi,--- , y m ) provided that one has the right structure on G (and therefore on F), namely 
for G 

(H3') The function G is locally Lipschitz continuous w.r.t. y = (y\, ■ ■ ■ ,y m ), t, p y and 
there exists I = (li, • ■ ■ , l m ) £ W 71 and C3, C4, C5 > such that, for all % — 1, • • • , m and 
for almost every £ = (x,y,t,p x ,p y ) £ K 2 ( n + m )+ 1 

\D yt G(0\ < Cz\p yt + k\ , \D t G(0\ < C 4 (l + \ Py \ + \G(0\) , \D Py G(0\ < C B . 

With this assumption, the proof of Theorem 11.11 remains the same : the only slight 
difference is that we have to do the estimates for all the D y .(w a (x,y,t) + kyi) and the 
first inequality for \D yi G(£)\ above implies that we can do it separately for each i. 

For the homogenization problems, we may have additional dependences in x, y, t and 
not only in the fast variables e~ l x, e^y, e~ l t. In this case, x, y, t become parameters in the 
ergodic-cell problem as are p x ,Py The main additional difficulties consists, on one hand, 
in obtaining the right dependence of F(x,y,t,p x ,p y ) in x,y,t in order to have a Strong 
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Comparison Result for the limiting equation and, on the other hand, in the convergence 
proof, to take care of the lack of uniform continuity of F in x, y, t. It is worth pointing 
out that this last difficulty arises since we do not have (a priori) Lipschitz continuous 
solutions for the ergodic-cell problem. 

For the first difficulty, since one can prove, in an analogous way as in the proof of 
Lemma EH1 that F(x,y,t,p x ,p y ) is continuous, then the dependence in x does not create 
any problem. But the dependence in t and y seems less clear and we do not see how to 
show that (typically) \D y F(x, y, t,p x ,p y )\ < (7(1 + \p y \) and an analogous property for the 
t-dependence. Here the fact not to have a Lipschitz continuous solution of the ergodic-cell 
problem creates a difficulty and we do not know how to solve it. 

On the contrary, the second difficulty can be solved by using a similar trick to the 
F fc -one and is not a real problem. 

Finally, even if, in this paper, we restrict ourselves to the case of first-order equations, 
there are easy extensions to second-order equations : for example, Theorem 11.11 extends 
without any difficulty if we add a u yy to the equation, which corresponds to homogeniza- 
tion problems with a vanishing viscosity term of the form £ 2 u £ yy . Such type of extensions 
will be considered in a future work. 



Appendix : Proof of (1131) 



Since (x, y, t) is a strict maximum point of u — <p, standard arguments show that, for e 
and S small enough, the function 

1 / -l / \ x ~ x I 2 \y ~ v I 2 1^ — ^ I 2 



U £ (x, y, t) - <f>(x, y, t) - eV(e- l x', e~W, e"Y) - 



6 2 S 2 6 2 



has a local maximum point near (x,y,t,x,y,t). We denote it by (x,y,t,x',y',t') for 
the sake of notational simplicity. Moreover we know that (x, y, t, x', y', t') converges to 
(x, y, t, x, y, t) as e, 5 — > and 

|<£ — x' I \y — y'\ 2 \t — t'P 

H — > as 5 — > for any e > 0. 



6 2 6 2 6 2 

We set 



2{x-x') 2(y-y') 2{t-f) 



2 



The viscosity inequalities for U e at (x, y, t) and for V at (e 1 x l ', e 1 y', e 1 t') (after applying 
the obvious unsealing) read 

4>t(x, y, t) + q t + F(e~ l x, e^y, e~h, D x <p(x, y, t) + q x , D y (p(x, y, t) + q y ) < , 

q t + F k (e' 1 x',e~ 1 y',£^H',q x + p x ,q y +p y ) > F k (p x ,p y ) 

where we recall that p x = D x <p(x,y,t) and p y = D y (f>(x,y,t). 

But, by the above recalled property on the penalization terms, we have 

\e^x — s~ l x'\ 2 + \e~H — e~Y| = s~ l 5o(l) where o(l) — > as e, 5 — > , 
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and, on the other hand, since (x,y,t) — > (x,y,t), we have 

| {D x (f)(x, y, t) + ft) - (ft. + p x ) | -> as e, S -> . 

Therefore if we choose 5 <C e and e, 5 small enough, we have 

\e- 1 x-e- 1 xf,\e-H-E-H'\,\(D x <f,(x,y,t) + q x )-(q x +p x )\ < AT 1 . (22) 

Next, we use this information together with (H3) to obtain 

F(e~ l x, e^y, e~H, D x (f)(x, y, t) + q x , D y (f)(x, y, t) + q y ) > 

F{e~ 1 x, e~ l y\ e~H, D x <f>{x, y, t) + q x , D y <j>(x, y, t) + q y ) - C 3 \e^y - £~V | \D y <f>(x, y, t)+q y \ 

-C 5 \D y (j)(x,y,t) - D y (/)(x,y,t)\ . 

And thanks again to the property on the penalization terms, we can choose 5 small enough 
in order to have 

\e~ y — s^y'l \D y <p(x, y,t) + q y \ = o(l) as e, 5 — > . 

Using this and ff22l yields the final estimate of F{e~ l x, e^y, e~H, D x <f)(x, y, t)+q x , D y (p(x, y, 
q y ), namely 

F k (e- 1 x', e~V> D x (f)(x, y, t) + q x , D y <f>(x, y, t) + q y ) - o(l) . 

In order to conclude, it suffices to subtract the above viscosity inequalities and to use 
this information : this yields 

4> t (x,y,t) - o(l) < —F (p x ,p y ) , 

and we let e, 5 tends to 0, choosing 5 sufficiently small compared to e. 

Remark 4.1 This detailed proof shows why the u F k -trick" can be useful : without using 
it, we would have to compare the terms F(e _1 x, e^y, e~H, D x (p(x, y, t) + q x , D y (p(x, y, t) + 
q y ) and F(e~ 1 x', e^y', e~ l t' , q x + p x ,q y + Vy)- As the above proof shows it, the "good" 
assumption (H3) allows to treat the differences in y and D y u and it seems to be a natural 
assumption to do it. But, for the x and t variables, we do not have such type of assumption 
and the u F k -trick" seems the only way to handle the corresponding differences. 

We have choosen to present this u F k -trick" in this way in order to point out the 
difference between having (H3) and not having it, but we could have defined F k in a 
different way in order that it takes care of all the variables. 
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